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We study the effect of surface adsorption of 27 different adatoms on the electronic and magnetic 
properties of monolayer black phosphorus using density functional theory. Choosing a few repre¬ 
sentative elements from each group, ranging from alkali metals (group I) to halogens (group VII), 
we calculate the band structure, density of states, magnetic moment and effective mass for the en¬ 
ergetically most stable location of the adatom on monolayer phosphorene. We predict that group I 
metals (Li, Na, K), and group III adatoms (Al, Ga, In) are effective in enhancing the n-type mobile 
carrier density, with group III adatoms resulting in lower effective mass of the electrons, and thus 
higher mobilities. Furthermore we find that the adatoms of transition metals Ti and Fe, produce a 
finite magnetic moment (1.87 and 2.31 hb) in monolayer phosphorene, with different band gap and 
electronic effective masses (and thus mobilities), which approximately differ by a factor of 10 for 
spin up and spin down electrons opening up the possibility for exploring spintronic applications. 


I. INTRODUCTION 

Thin body materials and devices are being considered 
as the most promising candidates for 5nm semiconduc¬ 
tor technology and beyond by semiconductor industry. 1 . 
This has lead us to actively study various 2D layered 
crystals such as graphene 2 , silicine 3,4 , germanene 5 , tran¬ 
sition metal dichalcogenides (M 0 S 2 , MoSe 2 , WSe 2 etc.) 6 , 
due to their excellent material properties 7-11 and re¬ 
duced short channel effects. Field effect transistors 
(FET’s) based on graphene 12 , and M 0 S 2 13 , have already 
been demonstrated with performance superior to con¬ 
ventional Si based FET’s. More recently, layered black 
phosphorous (dubbed as phosphorene) has joined this 
promising family of 2D crystals. In contrast to transi¬ 
tion metal dichalcogenides, the carrier (hole) mobility in 
phosphorene is quite high (^ 1000 cm 2 /V.s), which in 
turn gives a Ion/^-off current ratio of ~ 10 5 at room 
temperature 14-17 . 

Phosphorene is a layered material with individual lay¬ 
ers having a puckered structure and the layers being 
stacked together by a weak van der Waals interaction, 
allowing for mechanical exfoliation. In each layer, a phos¬ 
phorous atom covalently bonds with three other atoms, 
forming a honeycomb structure. Semiconducting phos¬ 
phorene has an intrinsic direct bandgap of ^ 0.31 — 1.0 
eV (bulk - monolayer), which is inversely proportional to 
the number of layers 15,18 . Because of it’s direct bandgap 
at T point, phosphorene is a potential candidate for ap¬ 
plications in the optoelectronic devices 14,15,19 21 . Most of 
the 2D materials have high electron mobility and display 
n-type device characteristics, while phosphorene also has 
a high hole mobility and FET’s having both p-type and n- 
type characteristics have been demonstrated 14,15 . Addi¬ 
tionally, phosphorene shows highly anisotropic behavior 
in effective mass and mobility, which can be controlled us¬ 
ing appropriate uniaxial or biaxial strain 22 . Black phos¬ 


phorus has also been proposed as an effective anode ma¬ 
terial for Li ion batteries 23 . Among other advantages, 
black phosphorous is known for it’s chemical inertness 
and superior transport properties 17,19,20,22,24-26 . 

For practical use of any semiconductor material in FET 
devices, doping is essential to enhance the carrier concen¬ 
tration and to control threshold voltage. Since 2D lay¬ 
ered materials have a very high surface area, adatom ad¬ 
sorption can be very effective strategy for carrier doping 
and band structure engineering. For example, adsorption 
of hydrogen/halogen, metal and molecules in graphene 
can change its electronic band structure, transforming 
graphene, from metal to semiconductor/insulator 27 29 . 
Adatom adsorption in 2D transition metal dichalco¬ 
genides creates donor or acceptor levels in between the 
valence and conduction band of the parent material, lead¬ 
ing to electron or hole doping of wide bandgap semicon¬ 
ductors like MoS 2 30 - 33 . Both experimental study and ab 
initio simulation of potassium adsorbtion in M 0 S 2 have 
shown to dope M 0 S 2 with n-type carriers 33,34 . Similar 
studies have been done for selected adatom adsorption 
in monolayer black phosphorous range of valences, in¬ 
cluding s and p valence metals, transition metals, and 
semiconductors, hydrogen and oxygen. 35,36 Further ox¬ 
idation of phosphorene has also been studied and may 
provide more stability and alternate route to synthesis of 
phosphorene. 37,38 . 

In this article, we use ab-initio spin polarized den¬ 
sity functional theory (DFT) calculations to systemati¬ 
cally explore the possibility of band structure engineering 
and doping possibilities of monolayer black phosphorous 
(MBP or monolayer phosphorene) via surface adsorption 
of various adatoms. We consider a few representative el¬ 
ements from several group in the periodic table, ranging 
from alkali metals (group I) to halogens (group VII) and 
our aim is to provide a comprehensive guideline for dop¬ 
ing strategies of MBP via surface adsorption of adatoms. 
In addition, we investigate the structural and magnetic 
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FIG. 1. (a) Top view and (b) side view of monolayer black phosphorous, (c) top view and (d) side view showing the five 
different adatom locations: light green at H (hexagon center), dark blue at T u (top of the phosphorous atom in the upper 
layer), dark green at Ti (top of phosphorous atom of lower layer), red at B u (bridge of the upper layer P-P bond) and finally 
dark brown at Bi (bridge of the bottom layer P-P bond), (e) Band structure of intrinsic MBP with dashed line showing the 
Fermi level, and inset displaying the symmetry points in the Brillouin zone, (f) Spin resolved density of states (DOS) for MBP. 
Red (t-spin) and blue (J^-spin) color lines are for complete DOS of MBP and black and purple color lines represent projected 
DOS of P atoms. 


properties, and the charge transfer from the adsorbed 
adatom to MBP and vice versa. Five different sites are 
chosen for adatom adsorption, as shown in 1(c) and 
(d). These are named as T u (on top of the phospho¬ 
rous atom in upper layer), B u (top of bridge of upper 
layer), H (centre of the hexagon), Ti (on top of phospho¬ 
rous atom in lower layer) and finally Bi (above the bridge 
in the lower layer). Our work broadly reaffirms the find¬ 
ings of Ref. [35] , and extends the study to include more 
adatoms such as group III elements (Al, Ga, In) among 
others, which are the most effective electron dopant with 
the lowest effective mass and possibly highest electron 
mobility. Additionally we also report the effective mass 
of spin polarized carriers, which is useful for exploring 
possible usage in electronic and spintronic applications. 


The paper is organized as follows: in Sec. II, we present 
the parameters and other details of the computational 
method used. In Sec. Ill, we briefly review the crystal 
structure and electronic properties of MBP. This is fol¬ 
lowed by systematic presentation of the results and dis¬ 
cussion in Sec. IV and finally we summarize our findings 
in Sec. V. 


II. METHODOLOGY AND SIMULATION 
DETAILS 

We use the spin polarized density functional the¬ 
ory (DFT) calculations, as implemented in the ATK 
package 39 ^ 41 for geometry relaxation, as well as the elec¬ 
tronic and magnetic property predictions. The exchange- 
correlation energy is treated by the generalized gradi¬ 
ent approximation (GGA) using the Perdew, Burke, and 
Ernzerhof (PBE) pseudo potential 42 . The Brillouin- 
zone integrations are performed with Monkhorst-Pack 
K-points grid with 8x8x1 Brillouin-zone sampling, along 
with energy cut off of 75 Hartree. We have used a double 
zeta polarized basis set for all the calculations performed 
in this work. We have used a 3x3x1 simple orthorhom¬ 
bic phosphorene super cell (containing 36 phosphorous 
atoms) with a vacuum region of 10 A along the c-axis, 
which is sufficient to eliminate interaction among mono- 
layers. 

The charge transfer from the adatom to the phospho¬ 
rene substrate is calculated using the Mulliken charge 
transfer method, i.e., by calculating the difference of va¬ 
lence electron charge of the isolated adatom and charge 
of the adatom in the vicinity of phosphorene; p* = 
Z a ~ Pa 43 - Here Z a indicates the valance electron charge 
of the isolated adatom with an appropriate sign, for eg. 
Mg has Z a = —2e where e ~ 1.6 x 10“ 19 eV, and p a 
is the calculated electron charge on the adatom over the 
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Additionally we also study the adatom mobility and 
diffusion kinetics on monolayer phosphorene using the 
nudge elastic band (NEB) approach which calculates the 
activation energy barriers encountered by the adatom in 
moving from one stable location to another. The NEB is 
fairly accurate approach to determine minimum energy 
path between initial and final positions 46,47 . For this, 
Broyden optimization scheme has been employed with 
eight intermediate images between initial and final posi¬ 
tions to find the potential energy barrier between them. 
We find that for all the twenty seven adatoms considered 
by us, the minimum potential energy barrier is 0.6 eV for 
moving from the most stable location to the next stable 
location. Since the minimum potential energy barrier of 
0.6 eV is much greater than the available thermal energy 
at room temperature (0.026 eV), this implies the stabil¬ 
ity of all considered adatom against diffusion from their 
stable location. 

Finally for further confirmation of the binding energy 
we redo the spin polarized virtual crystal relaxation cal¬ 
culations in Quantum Espresso package 44 taking van der 
Waals forces into account. This was necessitated on ac¬ 
count of the discrepancy in some of our ATK reported 
binding energies, with those reported in Refs. [35 and 
36]. For the binding energy verification in QE we use 
the projector augmented wave (PAW) pseudopotentials 
within the GGA-PBE approximation and the modified 
results of both the calculations are shown in Table I. 


III. PHOSPHORENE: CRYSTAL STRUCTURE 
AND ELECTRONIC PROPERTIES 


FIG. 2. Panels a), d), g) and j) show the band structure, 
panels b), e), h) and k) display the spin polarized total (red 
for t-spin and blue for N s pi n ) an d projected DOS (black for 
t-spin and magenta for .f-spin) of MBP with selected adatoms 
of group I, II, III, and IV, i.e. K, Mg, A1 and Si respectively. 
Panels c), f), i) and 1) show the top and side view of the 
optimized adatom location. 


monolayer phosphorene layer. Thus, p* < 0 (p* > 0) 
implies excess electron charge transferred from (to) the 
adatom to (from) the phosphorene monolayer. 

To establish the stability of each of the reported 
adatom location on phosphorene monolayer, we perform 
ab-initio molecular dynamics simulations using Quantum 
Espresso (QE) package 44 . The MBP-adatom supercell 
is equilibrated at 400 K temperature and the dynam¬ 
ics of the adatom is monitored for 1 ps. The kinetic 
energy cutoff of is taken to be 40 Ry. The exchange- 
correlation energy is treated using PBE pseudopotential 
and van der Waals correction is included under Grimme- 
D2 approximation 45 . Ionic temperature is controlled us¬ 
ing Anderson thermostat. We find that while all of the 
adatoms vibrate around their most stable locations (re¬ 
ported in table I), there is no significant shift or disloca¬ 
tion from the most stable location at 400 K. 


Black phosphorus (BP) is reported to be the most 
stable form among all the allotropes of phosphorus (for 
example red and white phosphorous among others) un¬ 
der normal conditions. BP has an orthorhombic crystal 
structure (space group: Cmca), having eight phospho¬ 
rous atoms in the unit cell 48 . Multiple layers are stacked 
together to form BP, although unlike graphite the con¬ 
stituent layers are not flat, but puckered in nature. Based 
on tight binding calculations BP was found to be a di¬ 
rect bandgap semiconductor with an energy gap of ~ 0.3 
eV at the Z point 49 , which was later confirmed via DFT 
calculations by Qiao et al 18 The authors also calculated 
the band structure of MBP and found it to be a direct 
bandgap semiconductor with an energy gap of ~ 1 eV at 
the T point. Our calculation yields a very similar result 
for pristine MBP [see 1 (e)-(f)] and this also validates the 
parameters chosen for simulation in this article, as de¬ 
scribed in Section II. 


IV. RESULTS AND DISCUSSIONS 

As shown in 1, we have considered five different adatom 
locations - H, T u , Tb, B t and Bi for the adatom. We start 
by initially placing all the adatoms studied at these five 
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FIG. 3. Panels a), d), g) and j) show the band structure, 
panels b), e), h) and k) display the spin polarized total (red 
for t-spin and blue for ^-spin) and projected DOS (black for 
t-spin and magenta for .f-spin) of MBP with selected metallic 
adatoms, i.e. Au, Pd, Fe and Ti respectively. Panels c), f), 
i) and 1) show the top and side view of the optimized adatom 
location. 


locations, and let the structure relax. The most stable 
configuration is determined based on maximum binding 
energy and is chosen for the purpose of presenting the 
electronic band-structure data. The binding energy is 
calculated according to the relation: 


-^binding — -E'MBP T -E'adatom -E'MBP+adatom 5 (1) 

where Fmbp is the total energy of the super cell of pris¬ 
tine MBP, -E/adatom is the energy of the isolated adatom 
in the ground state calculated for the same super cell 
with the same parameters, and ^MBP+adatom is the to¬ 
tal energy of the combined phosphorene-adatom system 
after structural relaxation. In I, we have reported the 
maximum binding energy value among the five different 
adatom locations (H, T u , Tp B u and Bi) and presented 
the corresponding electronic band structure, DOS and 
adatom projected density of states (PDOS) results in 2, 
3 and 4. Note that in these figures the dashed line dis¬ 
plays the Fermi level of the MBP-adatom system and the 



-2 h " -- 1 ' - 1 * 
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FIG. 4. Panels a), d), and g) show the band structure, panels 
b), e), and h) display the spin polarized total (red for T s pi n 
and blue for Jospin) and projected DOS (black for t-spin and 
magenta for .f-spin) of MBP with selected adatoms of group 
V, VI, and VII, i.e. As, O, and Cl respectively. Panels c), f), 
and i) show the top and side view of the optimized adatom 
location. 


band structure of MBP (black lines) has been shifted so 
as to align the conduction bands of pristine MBP and 
the MBP-adatom system. 

Let us begin by discussing the case of alkali and al¬ 
kaline earth metals. As shown in I, all of them prefer 
H site among the five possible locations and the vertical 
distance (d) from the upper phosphorous layer varies de¬ 
pending on the size of the adatom. Smaller atoms have 
higher binding energies, which implies stronger binding 
to the underlying MBP layer. In case of alkali metals 
(group I), adatom adsorption leads to n-type carrier dop¬ 
ing, because of which Fermi level shifts to the conduction 
band of MBP [see the case of K adsorption in 2(a)]. For 
each of the three group I adatoms studied here, MBP 
turns out to be electron doped; as expected because of 
highly electropositive nature of alkali metal atoms. Since 
the electronic states of K have negligible contribution to 
the total density of states (DOS) near the CBM [see 2(b)], 
the K adatom (and other alkali metals) preserves the 
original electronic structure of phosphorene, but shifts 
the Fermi level to the CBM by donating free mobile elec¬ 
trons, and hence enhances the electrical conductivity sig¬ 
nificantly. In case of group II, adatom adsorption yields 
significantly different effect. As shown in 2(d), after Mg 
adsorption a filled energy band is present in the band gap 
of MBP. Electronic states corresponding to this particu¬ 
lar band are localized near the adatom, which is clearly 
visible from the PDOS plot shown in 2(e). Thus, the 
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only way these electrons can take part in electronic con¬ 
duction is via excitation to the conduction band of MBP. 
However, since the band is located ^ 0.7 eV away from 
the CBM, group II adatoms (Mg, Ca studied here) are 
not going to be effective as electron donors. The PD OS 
plots [see 2(b) and (e)] reveal that the up and down spin 
states have mirror symmetry, which implies equal num¬ 
ber of up and down spin electrons and accordingly all the 
adsorbed atoms of group I and II are non-magnetic. 

Next, we consider the group III and IV adatoms. In 
case of group III, all the adatoms (Al, Ga, In) are located 
on the H site and they have similar binding energies [see 
I]. Among the group IV adatoms C prefers the H site, 
while Si and Ge are more stable on the Bi site and the 
binding energy is strongly dependent on the size of the 
adatom [see I] . Adatom-phosphorene distance d depends 
on the size of the corresponding adatom, as expected [see 
I] . In case of group III, as a result of adatom adsoprtion, 
the Fermi level has shifted to the conduction band of 
MBP, making it a n-type semiconductor [see 2(g)]. On 
the other hand adsorption of group IV adatoms leads to 
the formation of fully occupied midgap energy levels, lo¬ 
cated approximately in the middle of VBM and CBM [see 
2 (j)]. Thus, while group III adatoms are very effective for 
n-type carrier doping in MBP, group IV adatoms are not 
suitable at all for this purpose. We have not found mag¬ 
netic moment on any of the adatoms of group III and IV, 
which is corroborated by DOS plots of spin up and spin 
down states, as illustrated in 2(h) and (k). We emphasize 
here that only considering the charge transfer to deter¬ 
mine weather an adatom adsorption leads to electron or 
hole doping can be misleading, as in the case of group III 
elements [see I]. In case of group III adatoms, the changes 
in the band structure lead to a suppression of electronic 
states in the valance band, and an enhancement of DOS 
in the conduction band [see lh] . Also note that group III 
adatoms in MBP, leads to a lower effective mass of the 
electrons (^ 0.4m e ), which implies higher mobilities and 
consequently excellent transport properties. 

Now we consider the case of the transition metals. Ac¬ 
cording to the binding energy values reported in I, all 
of them get adsorbed preferably on the H site. Among 
them, only Au, Ti and Fe are found to be magnetic [see I] , 
which is corroborated by up and down spin DOS, shifted 
with respect to each other along the energy axis [see 3(b), 
(h) and (k)]. Interestingly, magnitude of the bandgap for 
up and down spin electrons are different, which can be 
exploited for spintronic applications. Also note that the 
effective electron mass for the up and down spin electrons 
in Fe and Ti, differ by a factor of 10 (approximately), 
which will also be reflected in their respective mobilities. 
Each one of the adatoms gives rise to the midgap energy 
levels [note that the Fermi energy lies in between VBM 
and CBM in 3(a), (d), (g) and (j)] and thus they are not 
effective either as n-type or as p-type dopant in MBP. 
We note however, that these trap states can speed up 
the carrier recombination process and can be useful for 
devices which require fast switch off times 50 . 


Finally we describe the effect of group V, VI and VII 
adatoms on electronic band structure of MBP. Based on 
the binding energy values shown in I, group V, VI and 
VII adatoms prefer Bp Ti and T u site, respectively. As 
expected the binding energy, as well as the value of d 
depends on the size of the adatom. As reported in I, 
smaller adatoms are closer to the MBP layer and bind 
more strongly than compared to the larger ones. In case 
of group V, adatom adsorption does not dope the un¬ 
derlying MBP layer with n-type carriers directly [see the 
case of As in 4(a)]. However, the energy gap (indirect) 
between the midgap state and the CBM is very small, 
such that these electrons can be excited to the conduc¬ 
tion band of MBP to increase the n-type carrier density. 
On the other hand, the bandgap remains unchanged by 
the adsorption of group VI adatoms [see the case of O in 
4(d)] and thus they are not effective for carrier doping in 
MBP. Among the group VII adatoms, while F is found 
to be non-magnetic, a small magnetic moment (~ 0.15 
fiB per adatom) is observed in case of Cl and Br. This 
is confirmed by the asymmetry of the DOS plots of up 
and down spin states in case of Cl adsorption, illustrated 
in 4(g). Note that, the bandgap of up and down spin 
electrons have different magnitude, which can be useful 
for spintronic device applications. Moreover for group 
VII adatoms, the Fermi energy lies close to the valance 
band, which may assist in making MBP hole-doped. 

Before concluding let us compare the reported binding 
(or adsorption) energy in previous works 35,36 and in this 
work. Kulish et al. in Ref. [35], did their calculations 
using QE, considering PAW method with PBE as the ex¬ 
change correlation functional, while ignoring dipole- as 
well as van der Waals- corrections. Hu et al. in Ref. [36], 
used VASP, with PAW-PBE functional, included dipole 
corrections and ignored van der Waals corrections. In 
this paper, we have used two different methods for cal¬ 
culating the binding energies: (1) using norm conserving 
pseudo potential with PBE functional in ATK, ignoring 
van der Waals forces and (2) using PAW-PBE functional 
in QE, taking van der Waals forces into account, respec¬ 
tively. We find that the binding energy depends signifi¬ 
cantly on type of pseudo potential used and weather or 
not the van der Waals correction is included. Further 
since our QE calculations use the same basis set and ex¬ 
change correlation functional as Ref. [35 and 36] , the QE 
calculated binding energies are also similar. 


V. CONCLUSION 

To summarize, in this article we have systematically 
explored the effect of surface adsorption of 27 adatoms 
from different groups of the periodic table, ranging from 
alkali metals (group I) to halogens (group VII), on elec¬ 
tronic, structural and magnetic properties of MBP. We 
find that among the adatoms considered, group I alkali 
metals and group III elements make MBP n-type, with 
the MBP-group III adatom system having lower elec- 
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TABLE I. Adatom, its preferred adsorption site [see 1 c) -d)], its binding energy (calculated using Eq. 1), vertical distance of 
the adatom d from the upper phosphorous layer a , magnetic moment per adatom (in the units of Bohr magneton, /jlb) in a 
3x3x1 supercell, electronic charge transfer and the electronic effective mass (in the units of electron rest mass, m e ) calculated 
at the r point in the T—X direction. 


Group 

Adatom 

Adsorption 

Site 

Binding energy b 
(eV) 

(i) , (ii) 

Height Magnetic moment 
d (nm) on adatom (hb) 

Charge 
transfer (e) 

Effective 
mass ( m e ) 

jst 

Li 

H 

1.922 , 2.545 

0.146 

0 

-0.024 

1.26 


Na 

H 

1.677 , 1.554 

0.200 

0 

-0.176 

1.24 


K 

H 

2.01 , 1.576 

0.258 

0 

-0.348 

1.27 

II nd 

Ca 

H 

2.193 , 2.513 

0.201 

0 

-0.056 

0.65 


Mg 

H 

0.944 , 1.088 

0.197 

0 

-0.170 

0.53 

III rd 

A1 

H 

2.370 , 3.186 

0.161 

0 

0.380 

0.37 


Ga 

H 

2.071 , 2.901 

0.172 

0 

0.448 

0.43 


In 

H 

1.970 , 3.124 

0.182 

0 

0.36 

0.43 

IV th 

C 

H 

5.578 , 6.330 

0.045 

0 

-0.01 

0.576 


Si 

B, 

3.130,4.284 

0.129 

0 

0.274 

0.337 


Ge 

Bi 

2.658 , 3.391 

0.148 

0 

0.14 

0.328 

V th 

N 

Bi 

3.270 , 4.822 

0.065 

0.056 

-0.028 

0.746(t) ,0.763(1) 


As 

Bi 

2.015 , 2.408 

0.143 

0 

0.036 

0.362 

VI th 

0 

Ti 

5.627 , 5.950 

0.092 

0 

0.136 

0.696 


s 

Ti 

3.368 , 3.894 

0.153 

0 

0.104 

0.650 

VII th 

F 

T u 

3.428 , 4.005 

0.181 

0 

0.206 

0.699(t),0.728(4.) 


Cl 

T u 

1.829 , 2.442 

0.218 

0.14 

0.124 

2.157(t), 1.881(4.) 


Br 

T u 

1.428 , 2.061 

0.231 

0.155 

0.149 

1.718(1),1.465(|) 

Nobel 

Cu 

H 

2.485 , 3.556 

0.113 

0 

0.214 

0.836 

Metals 

Ag 

H 

1.359 , 2.270 

0.146 

0 

0.046 

0.494 


Au 

H 

2.301 , 2.625 

0.150 

0.242 

0.064 

0.502(t),0.429(4) 

Transition 

Ti 

H 

3.858 , 4.519 

0.150 

1.874 

-0.056 

1.471(t),0.158(4) 


Fe 

H 

4.604 , 4.293 

0.112 

2.316 

-0.006 

2.825(4),0.324(4) 

metals 

Ni 

H 

4.330 , 5.524 

0.109 

0 

-0.042 

0.765 


Pd 

H 

3.567,4.547 

0.133 

0 

0.522 

0.602 


Pt 

H 

5.236 , 6.099 

0.121 

0 

0.208 

1.856 


Sn 

H 

2.252 , 2.932 

0.164 

0 

0.092 

0.305 


a Note that d is not the shortest distance from the nearest phosphorous atom 

b Calculated using: (i) Quantum Espresso package — PAW-PBE pseudopotential and vdW correction (ii) QuantumWise ATK — norm 
conseving basis set with PBE and no vdW correction 


tronic effective mass (^ 0.4m e ) which should result in 
better mobility and transport characteristics. Group II, 
group IV, and transition metals give rise to mid gap 
states, and are not effective as dopants. However, Au, 
Ti and Fe adatoms induce spin-polarization with mag¬ 
netic moments of 0.24, 1.87 and 2.31 fiB , respectively, in 
a 3 x 3 x 1 super cell. Moreover in the case of Fe and 
Ti adatoms the effective mass of the up and down elec¬ 
trons differ approximately by a factor of 10, are likely 
to have different mobilities and thus they may be very 
useful for spintronic applications. Amongst the adatoms 
of group V, VI, and VII, while none of them makes MBP 
either n-type or p-type directly, for the case of group 
VII elements the Fermi level lies close to VBM and it 
may be useful for increasing hole concentration. Finally 


we have also checked the stability of all the adatoms at 
their most stable location, against temperature fluctu¬ 
ations by doing a MD simulations, and against surface 
diffusion by calculating the activation energy barriers for 
adatoms moving to other locations. 
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